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ABSTRACT 


The refined mathematical model of the force created by the light 
pressure of the Sun has been used to compute the solar radiation pressure 
force acting on the Mariner 9 (Mariner Mars 1971) spacecraft, taking into 
account the reflectivity characteristics of all its components. The results 
have been compared with values obtained from Mariner 9 observations during 
the cruise phase and found to be in agreement within 0. 1% of the values. 


vi 
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I. INTRODUCTION 


In this study of the solar radiation force acting on the Mariner 9 
(Mariner Mars 1971) Mars orbiter, we shall use the refined mathematical 
model derived in Ref. 1, and apply the formulae to each component of the 
spacecraft, the reflectivity characteristics of which being known from 
ground measurements. 

The Mariner 9 Mars orbiter is a complex spacecraft consisting of 
a number of components of various shapes asymmetrically arranged about 
its roll axis, pointing toward the Sun (Fig. 1). The orientation of the 
spacecraft is maintained by the Sun and Canopus sensors so that a non- 
inertial spacecraft-fixed frame of reference can be defined in the following 
manner. Taking the Sun- spacec raft direction as the positive direction of 
the z-axis and denoting r as the heliocentric position vector of the space- 
craft, we have (Refs. 2 and 3) 


k 


r_ 

r 


where r = | r|, and k is the unit vector along the z-axis. 

We now define the unit vectur u^, along the direction in which the fixed 
star (Canopus) is viewed from the spacecraft. (Because of the presumably 
infinite distance of the star, this direction is the same as the geocentric 
direction of the star. ) Then the expressions 


and 


j 


u c x k 

u c Xk 


i = j Xk 


define the unit vectors along the y- and x-axes, respectively (Fig. 2). 
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We obtain the pitch- yaw-roll system of coordinate axes (Fig. 1) by 
rotating the above described frame of reference about the x-axis by the 


angle K = 180° + <|>, where (Fig. 3) 

<t> = 32. 2° 

The unit vectors along the axes of this system are given by 


(pitch) = - i cos 4> - j sin <J> 

(X-axis) 

(yaw) = i sin 4> - j cos <f> 

(Y-axis) 

(roll) = k 

( Z - axi s ) 


or, explicitly, 


e, = 


— cos <\> 

u~ — 

C sin y 


— cos 4> cot Y 


- ( u c x -f) 


sin 4> 
sin y 


= 


— sin <p 

- \i~ — : — — 

C sin Y 


sin cf> cot y - X “r”) 


cos <J> 
sin y 


e 3 


L 

r 


where 

Y = arccos (u^ • k) ^ 0 

Denoting u r as the projection vector of the vector u^. in the XY-plane, 
and counting the cone angle from the roll axis (Z) and the clock angle from 
the vector u f in the XY-plane in the positive direction, we can list the 
positions of the main spacecraft components, relative to the XYZ frame 
of reference (Ref. 4), as follows: 

•p 

(1) Four solar panels, 213.9X90.2 cm each, flat, perpendicular 
to the Z - axis. 
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(Z) Main body of the spacecraft, containing propulsion module, fuel 
tanks, thermal blanket, etc., wrapped in a protective cover, 
thus forming an approximately spherically shaped body of an 
average radius 78. 5 cm. 

(3) High-gain antenna in the form of a paraboloidal dish of radius 
50. 8 cm and depth 18. 0 cm. In the deployed position, the cone 
angle of the antenna's axis of symmetry is -38.3° and its 
clock angle is 80. 1°. 

(4) Medium-gain antenna, clock angle of its axis of symmetry 
55.0°, cone angle 158.0°. 

(5) Low-gain antenna, parallel to the Z-axis, radius 14 cm, 
depth 6. 4 cm. 

(6) Maneuver engine, approximately paraboloidal, radius 11 cm, 
depth 24 cm, parallel to the Z-axis. 

Other smaller components, illuminated by solar rays, will be mentioned 
later. 

II. THE SOLAR RADIATION FORCE 

Using the same notations as in Ref. 1, we can write the vectorial 
expression for the solar radiation force in the form of a surface integral 
over the illuminated surface area S: 


F 


K(r ) / B (0 ) dS - 


C 2 K(r)u 



cos 0 dS 


In this expression, 0 is the angle of incidence, given by 

cos 0 = u • N 


( 1 ) 


( 2 ) 
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where N is the unit vector along the local normal to the illuminated surface 
area S. Hence, 

dS = NdS 

u is the unit vector along the spacecraft-Sun direction, i. e. , 


u = 


£ 

r 


Denote y as the portion of reflected photons, (3 as the portion of 
photons reflected specularly, and k as a constant which depends on tem- 
peratures and emissivities of the front and back sides of the reflecting 
surface (Ref, 5), such that 


k = 1 


for adiabatic surfaces; 


4 4 

e T - e T 
F F B B 

e T 4 + e T 4 
F F B B 


for structures which have different temperatures on the front and back sur- 
faces (e^, and e g are emissivities and and are temperatures of the 
front and back surfaces, respectively); 


k 



for structures where the front and back surfaces have the same tempera- 
ture, i. e. , for which TV = T^. 

h Jd 

In further calculations, we shall assume that T = T for all com- 

r B 

ponents of the Mariner 9 spacecraft. 
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Having defined y, (3, and k, we can now write the explicit form of 
the function B(9 ): 

B(0) = cos 9 + 2C^ cos 0 (3) 

and the constants C^, and C^, which appear in Eqs. (1) and (3), are 

Cj = | [vd - P) + Ml - Y)] | 

c 2 = i - py > < 4 ) 

C 3 = p Y = 1 - c 2 I 

The function K(r), where r = |7| is the heliocentric distance of the 
spacecraft, is defined by 


K q R 

K(r) = (5) 

r 

where K c _ is the solar radiation constant: 
oiv 

K OT5 = 1 . 010 X 10 17 kg-m/s 2 (6) 

The values of reflectivity characteristics p, Y> e p> e g> an< 3 surface areas S for 
Mariner 9 spacecraft components are listed in Table 1 (Ref. 6). The 
quantities k, C ^ , C^, C^, derived from these values, are also listed in 
Table 1. 
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III. THE FLAT -SURFACE COMPONENTS OF THE 


MARINER 9 SPACECRAFT 

For flat surfaces which, due to the orientation of the roll axis of 
the spacecraft, are perpendicular to the direction of the incoming solar 
radiation. 


N = u = - — , 0=0 

r 

Hence, from Eqs. (1) and (3), 

F = K(r)(Cj + C 2 + 2C 3 )S-^ 

and the normal force is 

F n = P n K(r) (7) 

where, from Eqs. (4), 

P n = 1 (2 + ( 3 ) + ^ (1 - Y ) ( 8 ) 

The values of P are given in Table 1. The total value of P for flat sur- 
n n 

faces (solar panels and small parts which can be approximated by flat 
surfaces) is 

P n = 9. 698 m 2 (9) 
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IV. THE HIGH-GAIN ANTENNA 


The high-gain antenna has the form of a paraboloidal dish. The 
axis of symmetry of the dish is inclined to the roll axis of the spacecraft 
by an angle of 38,3° (Fig. 4). The clock angle of the axis of symmetry 
of the antenna is 80. 1°, which means that the angle between the line OC 
and the positive direction of the X-axis is X = 47. 1°, as shown in Fig. 4. 

The formulae for the components of the solar radiation force for 
the parabolic antenna, which we take from Ref. 1, are given in an antenna- 
fixed coordinate frame of reference, defined in the following manner. The 
z-axis of the system is the axis of symmetry of the antenna, so that the 
unit vector along this axis is given by 

( cos X sin a 
sin X sin a 
- cos a 



or, in other words, 

z = X cos X sin a + Y sin X sin a - Z cos a (10) 

The unit vector along the x-axis of the system is then defined by 

u X e 3 __ _ 

~e - — : , cos a = u • e >0 

x s in a z 


or 


e 


x 


( sin X 
- cos X 
0 
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and, finally, the unit vector along the y-axis is given by 


e = e x e 
y z x 


or 


e 

y 



Hence, the transformation between the components of the solar radiation 
force in the XYZ- system and the components in the xyz -system is given 
by 


F 


F 

F 



( sin X 
- cos X 
° 


- cos X cos a 

- sinX cos ^ 

- sin cx 



(ID 


or 



( 0. 7420 
-0. 6704 
° 


-0.5261 
-0. 5823 
- 0. 6198 



Following Ref. 1, we shall denote 6 as the diameter of the antenna 
reflector and f, as the depth of the dish and introduce the auxiliary angle ft 
defined by 


tan ft 


2 


£ 

6 
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The components of the solar radiation force are then 


F =0 
x 


F = K(r)(C 3 I 21 + Cjl 22 - Casino.) 


( 13 ) 


where 


F = - K(r)(C-I,. + C,I,, + C,I cos a) 

Z J J 1 i. L* C* 


2 2 

1^1 = -ir6 (1 + 2 cot ft In cos £2) sin 2 a 


22 


3 1 


32 


1 r 2 1 - cos ft ^ x 

3 1:6 1 + coT ft (2 + secn) Sma 


j -rr6 I 2 ^l-2cot 2 ftlncosft 

- (1 + 6 cot 2 £2 In cos £2) cos 2 nj 

0 ,2 cos £2 

2tt6 - — — cos a 

1 + cos Q 


\ 


(14) 


I = 2tt6 cos a 


For the Mariner 9 high-gain antenna, tan ft = 


0. 7087, tt 6 2 = 0. 8107 m 2 . 


and 


= - 0. 1496 m 


I 22 = - 0. 0548 m 


I 31 = 0. 8682 m 


I 32 = 0. 5717 m 


I cos n = 0. 4993 m 


I sin a = 0. 3943 m 
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F rom T able 1 , we find that C ^ - 0.47, C ^ — 0.80, - 0.20, and, 

therefore, the components of the force are 

F = 0 
x 

F = - 0. 371 K(r ) 

y 

F = - 0. 842 K(r) 
z 

Finally, from the transformation Eqs. (12), we obtain 

F x = - 0. 155 K(r) 

F y = 0. 171 K(r) 

F z = 0. 891 K(r) 

By analogy with Eq. (7), we can write that 

P v = - 0. 155 m 2 

P y = 0. 171 m 2 

P, = P =0. 891 m 2 
Z n 

V. THE MEDIUM-GAIN ANTENNA 

2 

For the medium-gain antenna, Table 1 yields = 0. 086 m , but no 
reflectivity characteristics are given. Consequently, we have to consider 
the antenna as a flat surface. The antenna is shaped in the form of a cir- 
cular cylinder of radius 4. 8 cm with a dish at the end of radius 12. 5 cm. 
The depth of the dish is approximately one-half of the protruding part of 
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the cylinder. The total length of the antenna is h = 0. 61 m. The surface 
area of the antenna is thus 


A p = 0. 120 m 2 


Due to the lack of data, only the normal component of the solar 
radiation force can be determined. The other two components, F^r and 
Fy> although different from zero, cannot be determined. The force F^ 
is, therefore, 

F z = F N = 0.086 K(r) (17) 

VI. THE LOW- GAIN ANTENNA 


The axis of symmetry of the low- gain antenna is parallel to the roll 
axis (Z-axis) of the spacecraft. Using the same denotations as in Eqs. (13) 
and (14), we find for the antenna reflector: 


a - 0 

tan = 0.914 


tt6 


2 


0. 062 m 


2 


I 

I 

I 


21 

31 

32 
I 


0. 090 m 2 
0. 053 m 2 
0. 062 m 2 
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Hence, P 

n 


0. 093 m , and 


F x = F Y = ° 

F N = F z = 0.093 K(r) 


( 18 ) 


VII. THE PROPULSION MODULE 

The surface area of the propulsion module is approximately a sphere 
of radius 0. 785 m. Since the top part of the module is shaded by the 
maneuver engine, the illuminated surface area of the module is a spherical 
belt between the latitudes 0 and 82° (Fig. 5). The equation of the surface is 

$>(X, Y, Z ) = X 2 + Y 2 + Z 2 - 6 2 = 0 

where 6 is the radius of the sphere. The unit vector along the local nor- 
mal on the surface is given by 


N = § rad$ = JL 

| grad $| 6 


where R =6. Hence, 


cos 0 = u • N 



Z 

6 


Introducing polar coordinates, 

X = p cos <f> 

Y = p sin 4 > 

dXdY = pdpd<f> 
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we have 


dXdY c d pd<f> 

dS = — = op 


N * k J7 Z - p 2 


and 


cos 0 - - 


j? - p 2 


z - - 1777 


Substituting these values into the expression for the solar radiation force, 
given by Eq. (1), we find, after the integration. 


F X = F Y = ° 


2 ‘ 1 yfb" - a" + C„ + C, 


F z = ir ( 6 - a )K(r) 


2 2 


36 


2 3 


-T* c 3) 


(19) 


where a is the radius of the maneuver engine at the top. Denoting 


sin p = — = 0. 140 
0 


we find, because = 1 , for the normal force 


2 / 2 ^*i 2 \ 2 

F z = F n = tt 6 K(r ) 1 1 + — — cos p - sin p } cos p 


) 
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or 


F z = f n = z ‘ 358 K(r) ) 

P = 2. 358 m z ) 

n / 

and, also, for the surface area illuminated by sunrays, 

Ap = 1. 936 m^ (with maneuver engine) 
VIII. THE MANEUVER ENGINE 


( 20 ) 


( 21 ) 


The axis of symmetry of the maneuver engine is parallel to the roll 
axis of the spacecraft. We assume that the engine nozzle has the shape of 
a paraboloid of revolution and once again apply the formulae (13) and (14). 
Here 


tan Q = 4. 364, & = 0 


Hence, we obtain 


and, accordingly, 


I 

I 

I 


21 

31 

32 


0. 012 m z 

0. 014 m Z 


F x f y - 0 


F z = F N = 0. 032 K(r ) 


P = 0. 032 m 
n 


( 22 ) 
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IX. THE TOTAL SOLAR RADIATION FORCE 


Summing up all the components of the solar radiation force acting 
on all parts of the spacecraft, given by Eqs. (9), (15), (17), (18), (20), 
and (22), we find 

F x = - 0. 155 K(r) 

F y = 0. 171 K(r) 

F z = F N = 13. 158 K(r ) 

I 2 2 Z~ 

and the total force Fg R = I F^ + F y + F^ is 

F gR = 13. 160 K(r) (24) 

At the average distance of Mars, this force is 

Fg R = 2. 558 X 10 5 newtons 

During the cruise phase, the weight of the spacecraft was W = 558. 8 
kg. Hence, its acceleration due to the solar radiation was 

a _ = 4.58X10 ^ km / s ^ 

X. COMPARISON WITH OBSERVATIONS 

The simplicistic solar radiation force model, currently used in the 
Jet Propulsion Laboratory's double precision orbit determination program, 
is given by 

F = (Gk + Gl+G T )A_K(r ) 

' r x y J ' P ' ' 
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where A ^ is the effective cross-sectional area of the spacecraft, and G , 

Jr 3C 

G , G are solve-for parameters that are to be determined from 

y r 

observations. 

Since G and G are small, their determination from observations 
x y 

is rather inaccurate. However, the value of G^ has been evaluated during 
the cruise phase of the spacecraft. Its value is 

G = 1.2275 (25) 

r 

A comparison of the two models yields 

P (TOTAL) 

G — A- 
r " Ap(TOTAL) 

We have already found in Eq. (23) that the total value of P n is 

P = 13. 158 m 2 
n 

We also find that the total surface area of the illuminated part of the space- 
craft is 

A p = 10. 708 m 2 

Hence, we find 

G = 1.2288 (26) 

r 

Comparing the observed value (25) with the computed value (26), we see 
that they agree within 0. 1% of the value of G^. 
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Table 1. Reflectivity characteristics of spacecraft components 




LOW GAIN ANTENNA 



Mariner 9 spacecraft 






Fig. 4. The high-gain antenna 
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TO SUN 



Fig. 5. Propulsion module 
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